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Abstract. We study in detail the influence of the nuclear form factor both on the Born cross section 
and on the Coulomb corrections to the photo-production of muon pairs off heavy nuclei (-yZ — » fi~ Z) 
and in heavy-ion collisions {ZZ — > ZZ/j, + fi~). Our findings indicate a number of issues which have not 
been sufficiently described as yet in the literature: (i) the use of a realistic form factor, based on the 
Fermi charge distribution for the nucleus, is absolutely indispensable for reliable theoretical predictions; 
(ii) we checked quantitatively that the equivalent photon approximation has a very good accuracy for 
the discussed processes; and (Hi) we present a leading logarithmic calculation of the Coulomb corrections 
which correspond to multi-photon exchange of the produced (jl with the nuclei. These corrections are 
found to be small (on the percent level). Our result justifies using the Born approximation for numerical 
simulations of the discussed process at the RHIC and LHC colliders. Finally, we calculate the total cross 
section for muon pair production at RHIC and LHC. 



1 Introduction 

Lepton pair production in ultra-relativistic nuclear colli- 
sions was discussed in numerous papers (see [TJ[SJ[3] for 
a review and references therein). For dcfinitcncss, wc re- 
strict ourselves to equal charge numbers of the nuclei 
Z\ = Z2 = Z and symmetric Lorentz factors 71 =72 = 7, 
for the RHIC and the LHC colliders with parameters given 
in Table ED 

In the present paper, we primarily consider the pro- 
duction of a muon pair, but for completeness and com- 
parison, we first recall some results for electron-positron 
(e + e~) pair production and therefore make a slight de- 
tour. The production of a single e + e~ pair in the Born 
approximation is described by the Feynman diagram of 
Fig-Hi t ne corresponding cross section was obtained many 
years ago [J]. Since the Born cross section is huge 

(see Table the e + e~ pair production can be a serious 
background for many experiments. It is also an impor- 
tant issue for the beam lifetime and luminosity of these 
colliders [5|. This means that various corrections to the 
Born cross section, as well as the cross section for 71-pair 
production, are of great importance. The subject is in- 
herently difficult; a number of controversial and incorrect 
statements in the literature have been clarified in Refs. [TJ 



Table 1. Cross sections for the production of light lepton pairs 
at modern colliders 



Collider 


Z 


7 <?Bo?n [kb] 


*££ [b] 


RHIC, Au-Au 


79 


108 36.0 


0.209 


LHC, Pb-Pb 


82 


3000 227 


2.46 



Since the parameter Za is not small (Za rs 0.6 for 
Au-Au and Pb-Pb collisions), the whole series in Za has 
to be summed in order to obtain the cross section with 
sufficient accuracy unless higher-order corrections are oth- 
erwise parametrically suppressed. The exact cross section 
for single pair production g\ can be represented as the 
sum of the Born value, the Coulomb correction, and of 
the unitarity correction, 



Cl = "Bo 



CCoul + Ounit • 



(1) 
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The Coulomb correction acoui corresponds to multi- 
photon exchange of the produced e ± with the nuclei 
(Fig. [2]); it was calculated in Ref. [BlllOj. The unitarity 
correction <j uni t corresponds to the exchange of light-by- 
light blocks between the nuclei (Fig. [3]); it was calculated 
It was found that the Coulomb corrections are 
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Fig. 1. Feynman diagram for the lepton pair production 
ZZ — > ZZl + l~'m the Born approximation (/ = e, /i) 
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Fig. 2. Feynman digram for the Coulomb correction 

about 10 % while the unitarity corrections are about two 
times smaller (see Table [2]) . In the last column of Table [2] 
is shown the result of Baltz (see Ref. [TT|) obtained by 
numerical calculations using formula for the cross section 
resulting from "exact solution of the semiclassical Dirac 
equations." In fact, the employed formulas allow to cal- 
culate the Coulomb correction in the leading logarithmic 
approximation only and this may account for the discrep- 
ancies of the results for RHIC indicated in Table [2] For 
the case of electron-positron pairs, the leading logarithmic 
approximation is insufficient because of the large absolute 
magnitude of the correction. 




r 



Fig. 3. Feynman diagram for the unitarity correction 



In this paper, we present detailed calculations related 
to muon pair production. This process may be easier to ob- 
serve experimentally than e + e~ pair production. It should 
be stressed that the calculational scheme, as well as, the 
final results for the /i + fi~ pair production are quite differ- 
ent from those for the e + e~ pair production. 

The principal issues related to muon pair production, 
including the problem of unitarity corrections, have been 
considered in Refs. [§1[T2"] . In particular, using simple esti- 
mates, it was pointed out that: (i) the Born contribution 



can be easily calculated using the equivalent photon ap- 
proximation (EPA) which has in our particular case a good 
accuracy; (ii) contrary to the e + e~ case, the Coulomb cor- 
rection is small in the muon case (on the level of a percent). 
The last statement is of principal importance because it 
justifies the validity of the Born approximation for event 
generators of this process at the RHIC and LHC colliders. 

In a recent paper [13], the conclusion (i) has been con- 
firmed, but the point (ii) has been questioned. Namely, 
in Ref. [13], it was found out that the Coulomb correc- 
tions to muon pair production are rather large: —22 % for 
RHIC and -17% for LHC. These results have been ob- 
tained using the same formulas as for the e + e~ case with 
the minor changes. Below, we present a new calculation of 
the Coulomb corrections for muon pair production in the 
leading logarithmic approximation (LLA) ; our result is in 
agreement with a previous numerically small estimate of 
the Coulomb corrections as given in Ref. [12], but it is in 
strong disagreement with the result of Ref. [13] . 

We would like to note that the above features (i) and 
(ii) are directly related to the fact that both the electro- 
magnetic form factors of the nuclei F(K 2 ), F(Q 2 ) and the 
cross section for the virtual block J*(k) + J*(q) — > /J + fi~ 
drop quickly with increasing photon virtualities K 2 = 
— k 2 > and Q 2 = —q 2 « q 2 > 0. However, the scale 
of this decrease is much less for the nuclear form factor 
than for the virtual 7* + 7* — > block (by 7*, we 

here denote a virtual as opposed to a real photon). 

As a rule, the calculation of muon pair production for 
nuclear collisions is very laborious (for example, the exact 
expression for the Born cross section even for the case of 
simplified form factors is an eight-fold integral). There- 
fore, it is convenient to check the main points of various 
approximations using the simpler process of muon photo- 
production. In this case we have a possibility to perform 
relatively easily both the exact and approximate calcula- 
tions and compare them. 

This paper is organized as follows. In Sec. [H we study 
in detail the photo-production of muon pair off heavy nu- 
clei jZ — > An exact calculation of the total Born 
cross section for arbitrary photon energy starting from 
the threshold is carried out. The use of a realistic form 
factor instead of simplified form F(Q 2 ) = 1/(1 + Q 2 / A 2 ) 
turns out to be critically important for moderate pho- 
ton energies. In Sec. [3] the validity of the EPA is stud- 
ied both for a realistic and for a simplified representation 
of the nuclear form factor. Coulomb corrections to photo- 
production of the muon pair are studied in Secs.lU Predic- 
tions for the RHIC and LHC colliders are given in Sec. [5] 
and we conclude with a summary in Sec. [6] Throughout 
the paper, we use a system of units in which c — 1, h — 1, 
a = e 2 j (he) ~ 1/137 and denote the muon and nuclear 
mass as m and M, respectively. 
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Table 2. 


Coulomb and unitarity corrections to the e e 


pair production 


Collider 


gg£ (Refs. EEH1) 


§^ (Refs. 021) 


( Ref - m 


RHIC, Au-Au 
LHC, Pb-Pb 


-10% 
-9.4% 


-5.0% 
-4.0% 


-17% 
-11% 



2 Form factor and Born-level pair 
photo-production 

2.1 Form factors and nuclear charge distributions 

We first recall basic formulas related to the realistic and 
simplified form factor representations for the colliding 
heavy nuclei which are central to our investigation. For 
the realistic form factor, we employ a Fermi-type nuclear 
charge distribution in the form (see Refs. [M] ) 



p(r) 



Pa 



1 + cxp [(r - R)/a] 



(2) 



with a = 2.30/(41n3) fm, R = 6.55 fm for Au (mass 
number A = 197) and R = 6.647 fm for Pb (A = 208). 
This leads to the mean squared radius 



/ira\ : 
V~R ) 



R. 



(3) 



with VV) = 5.4338 fm for gold and y/ir 2 ) = 5.5041 fm 
for lead. The latter numbers are in very good agreement 
with the experimental values \J (?" 2 ) exp = 5.4358 fm for 
5.5010 fm for lead. 



gold and 



' cxp 



The nuclear form factor is defined as 



F(q 2 ) = — J p{r) e-' lq r d 3 r , = / p{r) d 3 r , (4) 

where q 2 « Q 2 and q is the three-vector part of the pho- 
ton four-momentum q. Its behavior is shown on Fig. [4] it 
is seen that for Q 2 > 1/R 2 , the form factor drops quickly 
with the growth of Q 2 . On the other hand, the cross sec- 
tion for the virtual block 7* + 7* — » p + pT drops quickly 
with the growth of Q 2 at Q 2 > W 2 = (k+q) 2 > (2m) 2 [see 
Eq. (p~9j) below]. For further consideration it is important 
that 

1/i? 2 « (30 MeV) 2 < W 2 . (5) 



For the simplified form factor, we use an approxima- 
tion of a monopole form factor corresponding to an expo- 
nentially decreasing charge distribution 



F(Q 2 ) = 



1 



? -/A 2 



(6) 



Its behavior is also shown in Fig. [5] This approximate form 
of the form factor is used, for example, in Refs. 15,12. 
1 1 3] and enables to perform some calculations analytically. 




10 



QR 



Fig. 4. Realistic (solid line) and simplified (dashed line for 
A = 80 MeV and dot-dashed line for A = 90 MeV) form factors 
vs. QR for Au 



For the concrete calculations reported in Refs. [12, 13], the 
value 

yl = 80MeV (7) 

is used for lead and gold. In the calculations below we also 
use this value unless otherwise stated. Another possibility 
is to use the connection of A with the mean squared radius 



2 ), in this case 



A = 




90 MeV 



(8) 



for lead and gold. 

Looking at three curves in Fig. |U one can come to 
the conclusion that the difference between the two choices 
of the A parameter should be negligible. We will show, 
however, that a transition from the realistic form factor 
to the simplified one with A = 80 MeV or A = 90 MeV 
results in a change of the total cross section for the muon 
pair production at the RHIC collider on the level of 10 % 
or 20 %, respectively. 



2.2 Realistic form factor: exact result for the Born 
cross section 

We start to discuss the role of the form factor on the basis 
of muon pair production by a real photon with the energy 
u off the nucleus with charge Ze and mass M: 

7 (fc) + Z(P) - + trip-) + Z(P') . (9) 
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In the Born approximation, this process is described by 
the Feynman diagram of Fig. and the corresponding 
diagram is also contained as a block diagram within the 
Feynman diagram for pair production off heavy nuclei. We 
assume that a real photon with 4-momentum k and a vir- 
tual photon with 4-momentum q = P — P' and virtuality 



Q 2 = -q 2 >0 



(10) 



collide with each other and produce a [i + fi pair with the 
invariant mass squared 



W 



(k + q) 2 = 2kq- 



We also use the notations 



= (k + P) 2 = M 2 + 2uM , <T = 



2„,3 



Z 2 a 



11) 



(12) 



where m is the muon mass. So, to measures the incoming 
photon energy in the rest frame of the incoming nucleus. 




1 



~P+ 



,p> 



Fig. 5. Feynman diagram for the photo- production of muon 
pair in the Born approximation. The incoming virtual photon 
has momentum k, the invariant mass squared of the pair is 
W 2 = (k+q) 2 = (p + +p-) 2 . The four-momenta of the produced 
leptons are p± 



The exact cross section for muon pair production a^z 
can be split into the form 



C7Z = OBorn + CCoul , 



where <7Bom corresponds to the Born cross section, and 
the Coulomb correction ccoui corresponds to multi-photon 
exchange of the produced [P 1 with the nucleus (Fig.[6|). 



~P+ 



P* 



Fig. 6. Typical Feynman diagram for a higher-order Coulomb 
correction to the photo-production of a muon pair 



It is well known (see, for example, Ref. [TH]) that 
the exact Born cross section for the process ([TJ as well 
as for electro-production can be written in terms of two 
structure functions or two cross sections <7t{W 2 ,Q 2 ) and 
as(W 2 ,Q 2 ) for the virtual processes 77^ — > n + ^~ and 
77s — * respectively (here, 7 is a real initial pho- 

ton, while 7^, and 7^ denote the virtual transverse and 
scalar/longitudinal photons with helicity At = ±1 and 
As = 0, respectively): 



da 



Ll ,„. a =a T (W 2 ,Q 2 )dn T (W 2 ,Q 2 ) 

+ a s (W 2 ,Q 2 )dn s (W 2 ,Q 2 ). 



(14) 



The coefficients dnx and dn$ are called the number of 
transverse and scalar virtual photons (generated by the 
nucleus). The cross sections or and 175 can be found in 
Appendix E of the review [TB] : 



Ana 2 



W 2 + Q 2 



Am 2 W 2 - 8m 4 - 2Q 2 W 2 



(W 2 +Q 2 ) 2 



1 + 



Am 2 W 2 - AQ 2 W 2 



(W 2 + Q 2 



where 

v 



16na 2 Q 2 W 
(W 2 + Q 2 ) 3 



Am 2 



W 2 ' 



2m 2 
W 2 



L = 21n 



W 
2m 



(1 + v) 



(15) 
(16) 

(17) 



Let us note that 



as 



Ana 2 
~W 2 

16ira 2 Q 2 



W 4 



[1 + 0(Q 2 /W 2 )} , 
at Q 2 <c W 2 , 



(18) 



(I 3 ) and 



A-KOL 2 



I6vra 2 VF 2 
(Q 2 ) 2 



for Q 2 > W 2 . (19) 



The number of photons can be found in Sec. 6 and Ap- 
pendix D of Ref. [16] 



driT — 
dn s = 
where 



Z 2 a 

7T 

Z 2 a 



i-y- 



M 2 y 2 



2 2 dW 2 dQ 2 
w > w 2 -J- n2 n2 
dW 2 



y = 



kq W 2 + Q 2 



kp 



2u)M 



(21) 



Integrating the cross section (fT¥|) over Q 2 in the region 



Q 2 min <Q 2 < Qi 



where (see Problem 3 to § 68 in 
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Fig. 7. A Born cross section for the realistic (solid line) and Fig. 8. The function f(x) from Eq. (f29|) for the realistic form 
simplified (dashed line) form factors (photo-production on Au) factor (solid line) and f(x/{RA)) from Eq. (J36J) for the simpli- 
fied form factor (dashed line) (photo-production on Au) 



Ref. [T7]) 



O 2 ■ 

mm, max 



2B = 



B =p VB 2 - C , 

2M 6 - M 4 W 2 - M 2 W 4 
(2M 2 + W 2 )s 

2uj 2 -W 2 {1-uj/M) 



(22) 



s - 2M 2 - W 2 



1 + 2lo/M 



C = 



M 2 W 4 



(23) 



(24) 



and over W in the region 2m < W < to, we obtain the 
exact result for the Born cross section presented in Fig. [7J 



2.3 Simplified form factor: exact result for the Born 
cross section 

The exact result for the Born cross section for the case of 
a simplified form factor can be obtained using Eqs. (|14p . 
HU), (OH, (J20D with the form factor ©. The result is 
shown by the dashed line in Fig. [Jj It is seen that calcu- 
lations with the simplified form factor give an accuracy 
better than 10 %, 5 % and 2 % at oj > 3.5 GeV, 8 GeV 
and 50 GeV, respectively. At RHIC, the region near the 
"accuracy threshold" (2m < uj < 8 GeV) gives a numer- 
ically important contribution, which accounts for about 
10 -j- 20 % of the difference between cross sections with 
the realistic and simplified form factors. 



3 Approximations to Born-level pair 
photo-production 

3.1 Realistic form factor: equivalent photon 
approximation (EPA) 

Let us recall the usual schema of the EPA, but with the ad- 
dition of an accurate treatment of the nuclear form factor 



(see, for example, Ref. [16]). For the case of high-energy 
photons u> ^ 2m, the most important contribution to the 
photo-production cross section stems from photons with 
very small virtuality Q 2 -€i W 2 [we recall the definition 



of uo in Eq. (12]) and that Q 2 = -q 2 



It means that 



we can ignore the contribution of the scalar photons in 
Eq. (fl4]) and the dependence of <7t on Q 2 ; besides we can 
simplify the expression for dnx from Eq. (|20[) . As a result, 
we obtain the simple approximate (EPA) expression 



j EPA _ 
ua Born — °77 



(W 2 )dn^(W 2 ,Q 2 ), 



(25) 



where 



a^(W 2 ) = 



4-Ka 2 
~W 2 



4m 
~W 2 



8m 4 



Qi 



Z 2 a 

w 4 

Aw 2 ' 



4m 
W 2 



Q 



(26) 



2 

min 



Q 2 



2 2 dW 2 dQ 2 



F 2 (Q 



W 2 Q 2 



(27) 



The quantities v and L are defined in Eq. (|17[) . 

Integrating this spectrum over Q 2 , we obtain (the up- 
per limit of this integration can be set to be equal to infin- 
ity in a good approximation, due to the fast convergence 
of the integral at Q 2 > l/R 2 ): 



dn 7 (W 2 



Z?a_(W 2 R\ dW 2 



2lj 



W 2 



The function 



/;km£)? 



(28) 



(29) 
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1.2 



1.1 



< £ 

D. O 

cqcq 



1.0 



0.9 



0.8 



10 



w[GeV] 



n „A „SEPA /„. 
and <7 B orn /<7Born 



Fig. 9. Ratio CTBora/ cr Borii (solid line) 
(dashed and dash-dotted lines), where dBom is the total Born 
cross section from Eq. (| 14f) while OBom from Eq. (|32|l is cal- 
culated for the case of realistic form factor and erf from 
Eq. (|39[) is calculated for the case of a simplified form factor 
with A = 80 MeV (dashed line) and yl = 90 MeV (dot-dashed 
line) (photo-production on Au) 



is presented in Fig. [8l It is large for small values of . 



fix) = In 



Co 



for 



x < 1. 



(30) 



(The value of constant C depends slightly on the ratio 
a/R: we obtain Co = 0.166 for gold and Co = 0.163 for 
lead.) However, f(x) drops very quickly for large x, 



m < - 



for 



x > 1 



(31) 



Finally we obtain 
Z 2 a 



EPA 
"Born 



dW 2 ( W 2 R 
lm2 W* f { 2u 



a J7 (W 2 ) . (32) 



A comparison of this cross section with the exact result is 
shown in Fig. [HI It is seen that the EPA gives an accuracy 
better than 1 % already at uj > 1.3 GeV. This needs to be 
explained. 

Going from the exact expressions (|T3|). (fi~5]) , ([H?|). lgDJ) 
to the approximate ones l(23|) . (f2"6")l . (f2"5|) we omit terms 
of the relative order of Q 2 /W 2 , which are dropped be- 
fore the integration over Q 2 is done. After the integration 
with the "weight function" F 2 (Q 2 )/Q 2 the relative value 
of these corrections becomes of the order of 1/(R 2 W 2 ). 
In addition, the contribution of these correction terms is 
suppressed by a logarithmic factor. Indeed, the main con- 
tribution to the cross section in EPA is proportional to 
the large Weizsacker- Williams logarithm 



1/R 2 



dQ 2 

Q 2 



2 In 



2m 2 R J 



(33) 



while the omitted terms have no such a logarithm. There- 
fore, the actual parameter describing the suppression of 



the omitted terms to the differential cross section for pair 
production is numerically small indeed, 



1 



??EPA 



R 2 W 2 I W w 



(34) 



3.2 Simplified form factor: EPA 



The replacement of the realistic by the simplified form 
factor means that we have to replace the function / from 
Eq. (|29p by a function / which is obtained when we replace 
the form factor in the integrand in Eq. (|29p appropriately 
by the simplified nuclear form factor. The SEPA (S here 
stands for the simplified form factor) can thus be obtained 
using Eqs. (j2l)|) and (|2"6")) with the following expression for 
the number of equivalent photons, 



a m/2x Z 2 ^ r(W 2 \ dW 2 



(35) 



The function f[W 2 /(2cuA)] can be obtained analytically, 



f{x) = (l + 2x z ) ln U + 1 



2. 



(36) 



This is in contrast to f(x), which would be the equivalent 
of f(x) for a realistic form factor [see Eq. (|29p ]. Now, f(x) 
is large for small values of x, 



f{x) « In I A ) 2 



for 



x < 1 , 



but drops very quickly for large i: 



m < 



i 



6i 



7,4 



for 



x > 1. 



(37) 



(38) 



Its behavior is presented by the dashed line in Fig. [BJ 
where x — RA x. In view of the same leading logarithmic 
asymptotics for small argument [see Eqs. (l30l) and ([37)) ]. 
the functions / and / almost coincide for small values of 
x. 

Finally, we obtain for the simplified equivalent photon 
approximation (SEPA), 



SEPA _ 
"Born — 



Z l OL 



dW 2 ~Jw 2 \ ... 



Co JiwA/m 2 ) 



(39) 



For large photon energies, the function JiuiA/m 2 ) behaves 
as 



Tt \ 28 



i / \ 57 



for 



z > 1 



(40) 



A comparison of SEPA cross section with the exact result 
from Eq. (fT4]) is shown by the dashed and dot-dashed lines 
in Fig. [5] It is seen that the EPA gives a considerable 
better accuracy than the SEPA, once again confirming 
that the use of a realistic nuclear form factor is essential. 
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3.3 Realistic form factor: asymptotics for the Born 
cross section 



The high-energy asymptotic behavior of the Born cross 
section (for a realistic form factor) at large u> 3> 2m can 
easily be obtained using the EPA formulas ([25]). ([26]) . ([28]) 
with the asymptotic form of the function f(x) given in Eq. 
(j3TJj) . The final result is 



mula including the first correction ~ A 2 /(2m) 2 reads 



28 



a Born ~ If** 



In 



/ 2ui 



\ Km? 



43 
14 



1 



Co 



(41) 



It should be noted that the cross section (|4"Tj) provides a 
reasonable approximation only for large enough values of 
the photon energy to. Indeed, this cross section is positive 
only at 



1 



C = 4.4 GeV . (42) 



1 D 2 

to > w crit = - Km exp I — 



A comparison of the asymptotics with the exact Born cross 
section is given in Fig. [TU] It is seen that the accuracy of 
a simple expression (|4Tj) is better than 10 % only at very 
large u> > 20 GeV, showing that the realm of applicability 
of the high-energy asymptotics is limited. 




Fig. 10. Ratio a™ 
line) and ratio a™ Tn /aB 
Au) 



/o"Bom for the realistic form factor (solid 
(dashed line) (photo-production on 



3.4 Simplified form factor: result of Ivanov and 
Melnikov for asymptotics 



The cross section a 1 z in the high-energy limit was calcu- 
lated by Ivanov and Melnikov in Ref. [15] using the same 
expression © for the form factor of the nucleus and as- 
suming yl 2 /(2m) 2 <C 1. The corresponding analytical for- 



im _ TM 

a -yZ — °Born 



TM _ 28 

°Born — -g a 



IM 
7 Coul) 



28 



_im _ „ n 

°Coul - ~ "77" ^0 <^2, 



where 



12 f A 
35 V 2^ 



C* 2 = 0.928 {Zaf d 



(43) 



(44) 



(45) 



(46) 



We note that the parameter A 2 /(2m) 2 = 0.14 is small for 
muon pairs. A comparison of cr^ rll with the exact Born 
cross section ([T4")) is shown by dashed line in Fig. [TU1 

Two final remarks: (i) the SEPA asymptotics ([4"0)) is in 
accordance with the result of Ivanov and Melnikov ([44]) . 
as has already been noted in [T5]. (ii) The difference be- 
tween the high-energy asymptotics &^™ p for the realistic 
form factor (|4T]) as opposed to the high-energy asymp- 
totics cr™ ln for a simplified form factor is very small: 



IM 
"Born 



28 

0.012 —ao. 



(47) 



This is not surprising because the asymptotics are deter- 
mined by a region with small values of x — W 2 R/(2u>), in 
which the spectra of the equivalent photons for the real- 
istic and simplified form factors coincide (see Fig. [5} . 



4 Coulomb correction to the 
photo-production of pairs 

Having discussed the role of the nuclear form factor in the 
determination of the lepton pair production amplitude in 
the Born approximation, we now turn our attention to the 
role of Coulomb corrections. This is done according to our 
"master equation" ([13)) . The Coulomb correction is the 
leading correction beyond the Born amplitude, provided 
the latter is being evaluated with exact form factors. 

The Coulomb correction corresponds to Feynman di- 
agram of Fig. [6] The calculation of the Coulomb correc- 
tion for high photon energies (w 3> 2m) can be performed 
approximately using the result of Ivanov and Melnikov 
given in Eq. (|45]) . The ratio er™ ul /cBorn as presented at 
Fig [TT] is small. It is seen that the relative magnitude of 
the Coulomb correction is less than 1 % at to > 20 GeV. 
This is in accordance with the following estimate p~5]H2j . 

Due to the restriction of the transverse momenta of 
additionally exchanged photons to the range below A ~ 
1/R, the effective parameter of the perturbation series is 
not {Za} 2 but 

(Za f £ . (48) 

where W is the invariant mass of the muon pair. Besides, 
there is an additional logarithmic suppression because the 
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Fig. 11. Relative magnitude of the Coulomb correction 
(photo-production on Au) 



Coulomb corrections lack the large Weizsacker- Williams 
logarithm. Therefore, the actual parameter describing the 
relative value of the Coulomb correction is 



r?Coui = (Za) 2 



A 2 



W 2 Lww 



(49) 



which corresponds to Coulomb corrections of less then 1 % 
for oj > 20 GeV [we recall that Lww is defined in Eq. (f33|) ] . 
It is reassuring that the result of Ivanov and Melnikov 
confirms this estimate. Indeed the relative order of the 
Coulomb correction according to Eqs. (f44|) — ((46|) is 



IM 
7 Coul 



0.318 (Za) 



A_ 

2m 



\m 2 J 14 



?7Coul • 



(50) 



5 Predictions for the RHIC and LHC colliders 

We now turn our attention to the muon pair production 
in collisions of heavy nuclei. Let us therefore consider the 
process 

z(p 1 )+z(p 2 ) - ^( P+ )+ f r(p^+z(p{)+z(p^ . (5i) 

Its cross section can be calculated with a high accuracy 
by means of the EPA using the result (fT4|) for the exact 
cross section of the process, 



7 (fc) + Z(P 2 ) -» f i+( P+ ) + M "(p_) + Z{P' 2 ) . 



(52) 



For the RHIC collider, we use the parameters Z = 79 and 
7 = 100, the latter in order to be accordance with the 
value used in Ref. [13]. In the Born (B) approximation 
and with a realistic (Fermi, F) form factor, we have: 



zz 
"bf 



Z 2 a 



2 m 



JL 



^Bom^) = 0.193 barn. 

(53) 



In Eq. (|53|) . 7l = 27 is the Lorentz- factor of the first 
nucleus in the rest frame of the second nucleus; f(x) and 
crg orn (w) can be found in Eqs. (|2"9"j) and ([T3]). respectively. 
There is a 9.8% difference to the corresponding result for 
the simplified (S) form factor, still in the first Born ap- 
proximation, 

(54) 



erf I = 0.212 barn. 



This is in full agreement with the recent result 0.211 barn 
of Ref. [13J . The consistent use of A = 80 is crucial in 
order to obtain this agreement. We note in passing that 
A = 90 MeV results in a 22 % difference. 

A calculation for the Coulomb correction in LLA can 
be done using the result of Ivanov-Melnikov and taking 
into account Coulomb corrections to both nuclear lines 
(factor 2), 



zz 

7 Coul 



du) / ojR 

2m W V7L 



2(7 



IM 
Coul 



= -0.0072 barn. 



(55) 

It means that the relative value of the Coulomb correction 
is —3.7 % in full contrast to the recent result —22 % of 
Ref. [13], but in agreement with our parametric estimates. 

For the LHC collider, we use Z = 82, 7 = 2760, again 
in order to be in accordance with Ref. [13] . We have for a 
realistic form factor, 



zz 

'BF 



2.36 barn , 



and for a simplified form factor, 

off = 2.45 barn. 



(56) 



(57) 



This is in good agreement with the recent result 2.42 barn 
of Baltz [T3]. Again, an estimate for the Coulomb correc- 
tion can be obtained on the basis on an integration over 
the result of Ivanov and Melnikov, 



~zz 

°Coul 



-0.03 barn. 



(58) 



It means that the relative value of the Coulomb correction 
is —1.3 % in full contrast to the recent result —14 % of 
Baltz |13j . 

For completeness, we recall that in Table [U slightly 
different values were used for the relativistic Lorentz fac- 
tors at the modern colliders, namely, 7 = 108 (RHIC) and 
7 = 3000 (LHC) instead of 7 = 100 (RHIC) and 7 = 2760 
(LHC). In both cases, with the alternative values for 7 
and for a realistic nuclear form factor, we obtain results 
for <Tgp which are slightly larger than those in Eqs. ([53]) 
and (|5rJ|) . namely 0.209 barn for RHIC and 2.46 barn for 
the LHC (see Tabic p. 



6 Conclusions 

We have analyzed in detail the role of the nuclear form 
factor in the calculation of muon pair production cross 
sections in photon-nucleus and nucleus-nucleus collisions. 
At RHIC, the realistic (Fermi) nuclear charge distribution 
leads to predictions that deviate by 10 -j- 20% from the 
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corresponding values for simplified nuclear form factors. 
We also show quantitatively that the EPA is an excellent 
approximation to the muon photo-production for photon 
energies that exceed the rest mass of the produced pair 
(region lu 3> 2m) as well as for muon pair production at 
RHIC and LHC. 

We find that the Coulomb corrections for the muon 
production are less pronounced than for the e + e _ pair 
production. Our calculation in LLA leads to a decrease 
by about 1.3 -r- 3.7 % due to higher-order Coulomb effects 
at the LHC and RHIC colliders. 

Let us issue a few remarks regarding the obvious dis- 
crepancy of our results about the Coulomb corrections to 
those of the recent, interesting paper [13]. It is not obvious 
from the condensed presentation given in Ref. [13] whether 
or not the nuclear form factors have been taken into ac- 
count to all orders in Za. Therefore, the approach may 
need to be reexamined. Moreover, our parametric quanti- 
tative estimates given by the numerically small expansion 
parameter ([49| indicate that the Coulomb corrections for 
muon pair production should be smaller than those for 
e + e~ production. Coulomb corrections for the total pro- 
duction cross section of heavier lepton pairs would be even 
smaller, and in addition, we note that the Coulomb correc- 
tions also decrease with higher invariant mass W 2 . Under 
typical conditions, muons from the discussed process can 
be detected at the RHIC and LHC colliders with large 
values of W; this means that the Born approximation can 
be safely used in numerical simulations of this process. A 
correction on the order of 22 % for muons at RHIC, as 
obtained in Ref. [13], is larger than that for the e + e~ pro- 
duction and seems unrealistically large even if we allow 
for a large numerical prefactor multiplying the parameter 
([49]) . We also note that the calculation of the Coulomb 
corrections in Ref. [13] proceeds in the impact parameter 
representation. The numerical evaluation of integrals of 
this type is known to be notoriously problematic because 
of large numerical cancellations due to oscillations. In any 
case, a calculation of the Coulomb corrections beyond the 
leading logarithmic approximation is desirable. 

Finally, it should be mentioned than unitarity cor- 
rections to the muon production have been discussed in 
Ref. [12"1l9] with the following result: Unitarity corrections 
for the exclusive production of exactly one muon pair are 
large. However, the experimental study of the exclusive 
muon pair production seems to be a very difficult task, be- 
cause it requires that the muon pair should be registered 
without any electron-positron pair production, including 
e^- emitted at very small angles. The corresponding inclu- 
sive cross section is not affected by the unitarity correction 
and, indeed, close to the Born cross section. 
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